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The concept of magnetic connections is extended to non-ideal relativistic magnetohydrodynami- 
cal plasmas. Adopting a general set of equations for relativistic magnetohydrodynamics including 
thermal-inertial, thermal electromotive, Hall and current-inertia effects, we derive a new covariant 
connection equation showing the existence of generalized magnetofluid connections that are pre¬ 
served during the dissipationless plasma dynamics. These connections are intimately linked to a 
general antisymmetric tensor that unifies the electromagnetic and fluid fields, allowing the extension 
of the magnetic connection notion to a much broader concept. 
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In 1958 Newcomb El showed that in a plasma that 
satisfies the ideal Ohm’s law, two plasma elements con¬ 
nected by a magnetic field line at a given time will remain 
connected by a field line for all subsequent times. This 
occurs because the plasma moves with a transport veloc¬ 
ity which preserves the “magnetic connections” between 
plasma elements. This is one of the most fundamental 
and relevant ideas in plasma physics. 

The conservation of the magnetic connections imposes 
significant constraints on the plasma dynamics. More¬ 
over, it provides the very basis of central concepts in clas¬ 
sical (neither relativistic nor quantum) plasma physics, 
such as magnetic field line motion, magnetic topology 
and magnetic reconnection. However, as the plasma en¬ 
ergy increases (both in laboratory [3|-|5| and astrophysical 
contexts @-0), a generalization of these notions to the 
relativistic regime is required. Indeed, in a relativisic 
plasma different physical concepts must be taken into 
account, such as the distinction between magnetic and 
electric fields, which are reference-frame dependent. 

The key ideas to advance in the understanding of 
the magnetic connection notion in the relativistic regime 
have been first discussed by Newcomb [lj, and then 
thoughtfully investigated by Pegoraro [2]. In this latter 
work, using a covariant formulation for the dynamical 
equations of an ideal relativistic magnetohydrodynamic 
(MHD) plasma, it is shown that the magnetic connec¬ 
tions are preserved by taking advantage of the antisym¬ 
metry of the electromagnetic field tensor F^ v . In this 
case, the fundamental equation underlying the connec¬ 
tion concept takes the form p| 

J- ( dl= - (d"Up) ( dl a F a P) , (1) 

where d/dr = U v d v is the convective derivative along 
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the moving plasma with four-velocity XJ V and satisfying 
the ideal Ohm’s law U u F IJl ' = 0. Here, is the four- 
vector that connects two plasma elements. From Eq. 0 
it follows that if dl^F^" = 0 initially, it will remain null 
for all subsequent times. This generalizes the concept 
of magnetic connections, which can be recovered if one 
consider simultaneous events i'j . 

Despite these substantial advances in the understand¬ 
ing of the relativistic plasma dynamics, one should recall 
that there are plenty of non-ideal effects which play cru¬ 
cial roles in the dynamics of laboratory and astrophysical 
plasmas. In a non-ideal relativistic MHD plasma, effects 
such as thermal-inertial, Hall, or current inertia effects, 
modify the equation of motion and Ohm’s law. Therefore 
Eq. 0 is no longer valid and so the magnetic connec¬ 
tions are not conserved Hi,SHI- However, one might 
expect that in these dissipationless cases the magnetic 
connections could be replaced by more sofisticated “gen¬ 
eralized connections”. Along this work, we show that 
these generalized connections indeed exist. 

We suggest that the nature of the new connections can 
be understood through the emergence of a general anti¬ 
symmetric tensor field, which can be constructed as a 
combination of several physical quantities of the MHD 
plasma. This insight is deeply inspired by the ideas pro¬ 
posed by Bekenstein E2 and Mahajan E2, who have 
shown that the helicity invariants characterizing the ideal 
plasma dynamics can be generalized for relativistic plas¬ 
mas by realizing that fluid and electromagnetic fields cou¬ 
ple in one antisymmetric tensor. Similar ideas have been 
also successfully applied to develop relativistic theories 
of vorticity generation EMI, to study the topological 
constraints imposed by the plasma helicity E& E3 , and 
in the construction of generalized vorticity invariants in 
the framework of non-abelian plasmas |18j and spinning 
quantum plasmas Ell¬ 
in our analysis we intend to extend the connection con¬ 
cept by considering physical contents beyond the ideal 
MHD description. To this aim, we consider a plasma 
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governed by the generalized relativistic MHD equations 
recently derived by Koide 0. These equations re¬ 
tain many effects generally neglected in previous simpler 
models [2l], 0], such as thermal-inertial effects, thermal 
electromotive effects, current inertia effects and Hall ef¬ 
fects. The spacetime is assumed to be flat and defined 
by the Minkowski metric tensor r] /J/V = diag(—1,1,1,1). 
For an electron-ion plasma with density n , charge den¬ 
sity q = ne, normalized four-velocity t/ M (such that 
= —1) and normalized four-current density J 
the equations that govern the dynamics of the plasma 
are the continuity equation <9 M (gf7 M ) = 0, the generalized 
momentum equation 


law © in a form suitable to be antisymmetrized. We 
will see (below) that to achieve this, it is imperative to 
introduce an antisymmetric tensor that unifies electro¬ 
magnetic and fluid fields. 

First, using the continuity equation d l _ l (qU IJj ) = 0, 
and the charge conservation equation c/j = 0 (coming 
from Maxwell’s equations), we can rewrite the general¬ 
ized Ohm’s law ([3]) as 

j v iT + gujr J") - A nj v fr 

= n + qU v F> lv - A pJ v F» v + qW . (6) 


d v ( hU^U" + ) = -&*p + J V F^ , (2) 


(' 

and the generalized Ohm’s law 

d v 


—(tpj*' + j » u v ) - 

q q 2 


= + qU u F^ - A pJ v F» v + qRf 1 . (3) 


Here, h denotes the MHD enthalpy density, while n = 
pAp — A p, with p = p + +p- and A p = p + — p_ denoting 
the total pressure and the pressure difference between the 
fluids, respectively. We indicate with the subscript + (—) 
the positively (negatively) charged fluid. The reduced 
(rest) mass can be recognized as p = m+rri-/m 2 , with 
m = m+ + m_, wherease Ap = (m+ — m_)/m. The 
frictional four-force density between the fluids is 


R> i = -r l [Ji i +Q(l + Q)U li ], (4) 


Then, in analogy with the definition of the electromag¬ 
netic field tensor, we construct an antisymmetric flow- 
field tensor 

S IJV = d* (^U v ^j - d v , (7) 

and an antisymmetric current-field tensor 



The antisymmetric tensor 0 represents the covariant 
generalization of the fluid vorticity and it has been previ¬ 
ously introduced for relativistic hot plasmas in flat space- 
times m, 0 and in general relativity 0- In both 
above tensors, the four-vectors hU v /q and hJ v / q 2 may 
be viewed as the potentials (equivalent to the four-vector 
potential A M ) for the flow-field and the current-field, re¬ 
spectively. With these definitions Eq. 0 becomes 


where 0 is the thermal energy exchange rate from the 
negatively to the positively charged fluid, r/ is the plasma 
resistivity, and Q = U^J^. In Eq. ©, the current in¬ 
ertia effects arise from the left-hand side. Similarly, the 
thermal electromotive effects appear as inertial effects 
corrections in the left-hand side of Eq. ©^ whereas in 
the right-hand side it is taken into account the contri¬ 
butions of the thermal electromotive force and the Hall 
effect. 

As usual, F^ = d^A" — d v A 1 * is the electromagnetic 
field tensor ( A^ is the four-vector potential), which obeys 
Maxwell’s equations 

d v F^ = AnJ^ , 8 V F*^ = 0. (5) 


pJ^ + pqU v d» = 

+ qU„ {F^ + pA^) 

-J v {ApF^ - pS^ + pApA^) + qR^ . (9) 

To move forward in the antisymmetrization scheme, we 
define the general transport four-velocity 

W = U»- , ( 10 ) 

q 

which allow us to put the generalized Ohm’s law in the 
form 


Of course, F* IJV = (1/2 )e^ a0 F a p is the dual of F^, and 
indicates the Levi-Civita symbol. 

Although the inclusion of different relativistic effects 
complicates the MHD equations and a generalization of 
Newcomb’s results seems difficult to achieve, a hint may 
be gleaned directly from Pegoraro’s analysis Q, where 
the antisymmetric form of tensors is what allows to prove 
the existence of preserved magnetic connections. There¬ 
fore, the initial step of our search for generalized connec¬ 
tions consists in re-formulating the generalized Ohm’s 


pQ 


-cF 


-^[3 


+ ^Q = 


—d^n + u v (f^ + MH + -J v s^ + R» (11) 
2 q q 


Finally, we can manipulate the left-hand side of this equa¬ 
tion to reformulate the generalized Ohm’s law in the con¬ 
cise and elegant form 


£ M = U V M ,,V + R 1 ', 


( 12 ) 
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where we have defined the antisymmetric tensor 

MP" = - -¥-W^ , (13) 

A/x v ' 

with the antisymmetric general flow-field tensor 
W^ - A/xA'"' = d " ~ d v > ( 14 ) 

and E^ = d M [phQ/q 2 + ph/(qAp)] + (/x/A/x)x M , with 

= (is) 

V 9 / 9 \qj 2fiq 

The tensor Ad 111 ' defined in Eq. m is a generalized 
magnetofluid field tensor. This quantity represents an ef¬ 
fective field tensor which unifies the electromagnetic and 
fluid forces (through and W^ u ). Note that the ten¬ 
sor kW" contains the information of both the flow and 
current fields. We will show later that the effective field 
tensor A4^ v is crucial in revealing fundamental proper¬ 
ties of the system. This tensor is a more complex ver¬ 
sion of the unified magnetofluid field tensor introduced 
in Ref. [l2j for the simpler case of one-species hot rela¬ 
tivistic plasmas. 

To find the evolution of the generalized magnetofluid 
field tensor, we have to take the curl of Eq. m■ Apply¬ 
ing first e a| g 7A ,3 7 and then e a/3A< ^ to Eq. (fT2|) , we obtain 
the equation 

rj AA X(f) . . 

—— = d x U v M (l ’ v -d <l, U v M Xv —f-Z x ‘ l> +d x R' l ’-d <l ’R x , 
dr A/i 

(16) 

where now we have defined a general convective deriva¬ 
tive d/dr = U v d v , and introduced the antisymmetric ten¬ 
sor 


Z x<t> = qX^ _ q4 > x a _ ( 17 ) 

Therefore, Eq. (flTH) is an antisymmetric form of the gen¬ 
eralized Ohm’s law. 

Before proceeding further in showing the existence of 
generalized magnetofluid connections, let us analyze the 
physical content of the tensor field Z A A To this aim, we 
rewrite Eq. @ with the help of the continuity equation 
as 

u v d v (~uA = ~^J v d v (4 jA - -d*p + -J V F+ V , 

V 9 / 9 W J 9 9 

(18) 

and then we use it to cast Eq. m in a sum of different 
relativistic contributions to Ohm’s law 


_ 2^^ -\- 2^^ H - 2^^ 2^^ 


where 

T'Xfj) _ 
- 


Z = 


•7^0 _ 

£, u ~ 


Z M = 


A/x 


d x 




d x q 


d 4> 


A/x 

2/x 


n 





d * 


-JuF 


Xu 


d x 

i -j°d a (AjA 

+ d </> 

-J a d a (4 J A 


1 9 W )\ 


L 9 W )\ 


( 20 ) 


The antisymmetric tensors Z^f and Z x<t> are due to the 
thermal-inertial and thermal electromotive effects of the 
MHD plasma. In particular, Z x ^ contains the pressure 
contribution, which is absent if the equation of state of 
the MHD plasma is such that p , A p oc q. The contri¬ 
butions coming from the Hall effect in the generalized 

Ohm’s law are instead retained by the tensor Z X f , while 
Z M 

appears owing to current inertia effects. 

We are now able to demonstrate that when the fric¬ 
tional four-force density is negligible, e.g. when the 
evolution of the system is fast compared to the dissipa¬ 
tion time scale, there exist generalized magnetofluid con¬ 
nections that are preserved during the dynamics of the 
relativistic plasma. 

With this aim in mind, we define a general displace¬ 
ment four-vector Ax^ of a general element that is trans¬ 
ported by the general four-velocity 


Ax» 

At 


W + . 

A/x 


( 21 ) 


where At is the variation of the proper time and T)^ is 
a four-vector field which satisfies the equation 


M v<t, d x D v - M" x d^V v = Z^ . (22) 


The four-vector 2?^ contains all the (inertial-thermal- 
current-Hall) information of Z ^ v , and it turns out to be 
essential to prove the existence of the generalized con¬ 
nections. Similarly to Ref. @, we introduce the space¬ 
like event-separation four-vector dF = x ,Al — between 
two different elements. Two events are simultaneous in a 
frame where dl° = 0. From the definition of dF and the 
general four-velocity defined in Eq. m, it follows that 
this four-vector is transported according to ( d/dr)dF = 
U'v + (/x/A/x)©'" _ - (/x/A/x)^ = W{x a + dl a ) + 

(p/ Ap)V^(x a + dl a ) — U^{x a ) — (p/ Ap)V IJ, (x a ). There¬ 
fore, the four-vector dF fulfills 


4 dF = dl a d a (w + -r-vA . (23) 

ar \ A/x J 

With the help of Eq. (E3l) , and neglecting the frictional 
four-force density in Eq. m , we can calculate 

4(< H\M X *) = -d*U v (dl x M Xv ) 

+~T—dl\ (d x V u M 1 ' 4, - Z X(t> ) .(24) 
A (i 


(19) 
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This equation, which is in terms of the new four-vector 
field can finally be rewritten using Eq. m as 

^ (dhM^) = - (< ilxM xv ) 8* (u v + ■£j-V l ?j . (25) 

This is a crucial equation which reveals the existence 
of generalized magnetofluid connections that are pre¬ 
served during the plasma dynamics. Indeed, from this 
equation it follows that if dl\M. x< P = 0 initially, then 
d/dT(dl\M x< ^) = 0 for every time, and so dl\A4 x ^ will 
remain null at all times. Of course, regularity properties 
of the velocity field (1711) are assumed. 

The “magnetofluid connection equation” (1751) general¬ 
izes Eq. m for a relativistic electron-ion MHD plasma 
with thermal-inertial, Hall, thermal electromotive and 
current inertia effects. In particular, since 

dl x M M = dl x F** - -j^dlxW^ , (26) 

the generalized magnetofluid connections reduce to the 
well-known magnetic connections if (fi/Afi)dlxW x< ^ 0 
and dlo = 0. On the contrary, in the case F x ‘^ —>• 0, 
W x< ^ —> S x ^ and dlo = 0, the preserved connections are 
those of the fluid vorticity. 

Note that the four-vector connections (® are trans¬ 
ported by the general four-velocity U ' M + (n/A/j,)V^. The 
four-vector is fundamental to prove the existence of 
the generalized connections and it can be found by solv¬ 
ing the differential equation 

(27) 

obtained from Eq. (1771) . Here, A/^„ is the inverse of the 
matrix = S^v). The field exists 

since is non-sigular for non-ideal MHD. To show 

this, it is convenient to separate Eq. nil into time and 
spatial components. Defining the components of a gen¬ 
eralized electric-like field £ l = A4°* and a generalized 
magnetic-like field B k = (l/2)e ufc A4y, the determinant 

of is ||A4|| = ■ Bj . On the other hand, the 

time component of Eq. (1171) becomes U ■ £ = E° — R°, 
whereas the spatial components produce the equation 
U°£+UxB = E — R, implying that U°£-B = ^E — R^j-B. 

Therefore, as long as £ ■ B = E ■ B/U° — R ■ B/U a ^ 0, 
the tensor exists (i.e., ||Ad|| ^ 0). In the particular 
case of vanishing resistivity, it can be £ ■ B ^ 0 due to 
the non-ideal effects in E. 

We also observe that in general, if there is simultaneity 
bewteen two events, the condition dlxM Xcl> = 0 can be 
written as the vectorial conditions 

dl ■ £ = 0 and dl x B = 0 . (28) 

On the other hand, when the events are not simultaneous, 
Pegoraro has shown that simultaneity can be recovered 
resetting the time [2j. This can be achieved projecting 


the new trajectories of the fluid elements in a 3D space 
by changing dZ^ —> dl'^ = dl^ + ( U^ + fiH^/A[i)d\, such 
that in this new reference frame dl' 0 = 0. This allows 
to have unaltered generalized connections whenever H^ 
fulfills the generalized Ohm’s law (1171) . Thus, a possible 
solution can be expressed as H a = (A/i//Lt)7V’ Ma E A ‘ = 
A/)iaX M + Rffiad 1 ' ( AfihQ/q 2 + h/q). Notice the resemble 
with the V a field, which can be calculated from Eq. ED 
as 


2? a — NpotX^ + £ a ) (29) 

with ■ 

Our procedure has lead us to a completely general re¬ 
sult, and of course, also nonrelativistic systems are cov¬ 
ered. As an illustration, and for the sake of simplic¬ 
ity, let’s consider a non-relativistic Hall MHD plasma. 
For this case, we neglect the thermal-inertial and ther¬ 
mal electromotive effects. In the non-relativistic limit 
= 0{v 2 ) and Z X4> is negligible. Thus, the 
only contribution of the Z x< ^ field is Z^f, implying that 
X 4, = J^F^^/q ss because of our approxima¬ 
tions. Therefore, Eq. (1771) becomes (J v /q) + 

N’^F^Ja/q. A simple analytical solution of this 
equation may be obtained for quasi-uniform electromag¬ 
netic fields. Under this approximation, we have « 
— J v /q . Substituting this solution into Eq. (1211) . we can 
see that in this limit, the general vectorial velocity 

(30) 


preserves the topology of the non-relativistic limit (with 
no thermal-inertial effects) of the generalized magnetic- 
like field 


B ^ B — y—V X 
Aji 



(31) 


Therfore, in the considered limit the preserved connec¬ 
tions are related to the rotation of the fluid electron 
canonical momentum, as shown by Pegoraro [23]. We 
also note that it is common to assume m+ m_, which 
yields ^ ~ m-/m+ <C 1 and A/i « 1. In this case 
the transport velocity (1501) becomes the electron velocity, 
while the preserved connections are simply the magnetic 
connections. 

A difference of the theorvdeveloped here with the ideal 
MHD connection theorem [l],L3] is worth to be mentioned. 
In the ideal case, the magnetic connection concept is well 
defined when \E\ < \B\ (where U v is time-like). On the 
contrary, when \E\ > \B\, which may occur in extreme 
astrophysical environments [HI, [25|, the MHD approxi¬ 
mation breaks down since U v becomes space-like [lj. In 
the non-ideal case treated here, the same analysis is not 
directly applicable. The field strengths to be compared 
in the non-ideal system are |£| and |S| instead of \E\ and 
|J3|. In this case, the generalized Ohm’s law (1X71) gives 
that U = U°{£ x B)/\B\ 2 + (U ■ B/\B\ 2 )B + B x E, when 
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the frictional four-force density is neglected. Therefore, 
since the fields B and 8 depend on the general transport 
four-velocity, it is not straightforward to establish when 
the general transport velocity U is time-like or space-like. 
However, as long as U v is time-like, the generalized mag¬ 
netofluid connection concept is well defined. In this case 
it is also possible to have generalized magnetic-like field 
connections B in a frame where the magnetic field is null 
(.B = 0) because of the current and velocity fields. 

The generalization of the connection concept and of 
the connection equation to non-ideal relativistic MHD 
plasmas provides a strong theoretical framework for in¬ 
vestigating high-energy plasmas. In this extended frame¬ 
work, the generalized magnetofluid field tensor (fT 3 l) and 
the general four-velocity (1211) play the role that the elec¬ 
tromagnetic field tensor and the fluid four-velocity have 
in the ideal relativistic MHD description. Thereby, when 
dissipationless non-ideal effects are included, the pre¬ 
served connections are no longer related to the electro¬ 
magnetic field tensor F^ but to the magnetofluid field 
tensor A4 Ml '. This implies that it is possible to have re¬ 
connection of the magnetic field lines B while the gener¬ 
alized magnetic-like field lines B remain conserved (in a 


frame where dlo = 0). 

The conservation of these generalized magnetofluid 
connections profoundly affect the relativistic plasma dy¬ 
namics by forbidding transitions between configurations 
with different magnetofluid connectivity. Also, because 
of these constraints, they could be crucial to understand 
the formation of small scale structures resulting from a 
complex nonlinear dynamics. Therefore, it would not be 
surprising that the ideas presented in this letter will allow 
us to gain a more detailed comprehension of relativistic 
phenomena in plasmas. Furthermore, from the computa¬ 
tional point of view, these ideas can be applied to verify 
the accuracy of relativistic MHD numerical simulations, 
using the magnetofluid connection equation (12511 as guid¬ 
ance for reliability. 

Finally, the results presented here tell us that the con¬ 
nection idea, conceived first by Newcomb, could be a 
broader concept that permeates the entire energy scale 
in plasma physics. 
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